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ABSTRACT 

Context. Blue Compact Dwarf (BCD) galaxies constitute the ideal laboratories to test the interplay between massive star formation 
and the surrounding gas. As one of the nearest BCD galaxies, NGC 5253 was previously studied with the aim to elucidate in detail 
the starburst interaction processes. Some open issues regarding the properties of its ionized gas still remain to be addressed. 
Aims. The 2D structure of the main physical and chemical properties of the ionized gas in the core of NGC 5253 has been studied. 
Methods. Optical integral field spectroscopy (IFS) data has been obtained with FLAMES Argus and lower resolution gratings of the 
Giraffe spectrograph. 

Results. We derived 2D maps for different tracers of electron density (n e ), electron temperature (T e ) and ionization degree. The maps 
for n e as traced by [O n], [S n], [Fe m], and [Anv] line ratios are compatible with a 3D stratified view of the nebula with the highest 
n e in the innermost layers and a decrease of n e outwards. To our knowledge, this is the first time that a T e map based on [S n] lines 
for an extragalactic object is presented. The joint interpretation of the r e ([Sn]) and r e ([Om]) maps is consistent with a T e structure 
in 3D with higher temperatures close to the main ionizing source surrounded by a colder and more diffuse component. The highest 
ionization degree is found at the peak of emission for the gas with relatively high ionization in the main GH nR and lower ionization 
degree delineating the more extended diffuse component. We derived abundances of oxygen, neon, argon, and nitrogen. Abundances 
for O, Ne and Ar are constant over the mapped area within <0.1 dex. The mean 12 + log(0/H) is 8.26 ± 0.04 while the relative 
abundances of log(2V/0), log(Ne/0) and log(Ar/0) were ~ -1.32 ± 0.05, -0.65 ± 0.03 and -2.33 ± 0.06, respectively. There are two 

locations with enhanced N/O. The first (log(/V/0) 0.95) occupies an area of about 80 pcx35 pc and is associated to two super star 

clusters. The second (log(AVO) 1-17), reported here for the first time, is associated to two moderately massive (2 - 4 x 10 4 M Q ) 

and relatively old (~ 10 Myr) clusters. A comparison of the N/O map with those produced by strong line methods supports the use 
of N202 over N2S2 in the search for chemical inhomogeneities within a galaxy. The results on the localized nitrogen enhancement 
were used to compile and discuss the factors that affect the complex relationship between Wolf-Rayet stars and N / O excess. 

Key words. Galaxies: starburst — Galaxies: dwarf — Galaxies: individual, NGC 5253 — Galaxies: ISM — Galaxies: abundances 
— Galaxies: kinematics and dynamics 

1. Introduction ing c aused by a former encounter with M 83 (Ivan den Bergh 
, ™ ^ c , . mrrn ri980h . A wealth of studies at different wavelengths shows that 
Local Blue Compact Dwarf galaxies (BCDs, see ga i axy is peC uliar in several aspects. In particular, high reso- 
| Kunth&Osthn| | 2000t for a review) constitute ideal laboratories lution multiband photometry with the Hubble Space Telescope 
to test how the interaction between massive star formation gas (HST) revealed ^ existence of several re i at i ve ly young star 
and dust affects galaxy evolution Specifically, both winds of [q ^ central i(m of ^ galaxy wkh ^ masses of 
massive stars and supernova explosions i) inject mechanical ^ 2 _ 120 x 1Q 3 M ( Harril1jSD[2004). Among those, two 
energy to the interstellar medium (ISM), redistributing the gas massiye ^ : _ 2x 10 6 Mq) s Star Clusters (SSCs) at the 
q , within the galaxy and therefore, quenching (or igniting) future nudeus of ^ , separated by ^4 , Alonso-Herrero et al. 
_. star formation; 11) eject processed material into the ISM thus 2QQ4) and associated with a dense ct H n j onl 
causing a chemical enrichment of the galaxy detected ifl ^ radio at L3 cm and 2 cm (Xumer et al 2000)> 
NGC 5253 is an example BCD particularly suited for the st£md Qut ag candidates t0 be the t globular cluster(s) 
study of the interaction between gas and massive star formation, obseryed (Goriianetal 2001) Also? Lopez-Sanchez et al. 
since it is relatively close and a wealth of ancillaiy information is gjj showed how Ms al does not satisfy ±e Schmidt . 
available in all spectral ranges from X-ray to radio This galaxy Kennicutt law of star formation and seems to be slightly metal- 
belongs to the Centaurus A / M 83 complex ( |Karachentsev et alJ defident ^ arison with starbursts of ^ baryonic mass . 
\2MB> and is suffering a recent burst of star formation. The ex- Rut babl the most we ll. kn own peculiai'ity of NGC 5253 is 
i stence of an Hi plume exten ding along the optical minor axis ^ existence of areas with enhanced abundance of nitrogen. This 
( |Kobulnicky&Skillman|2 008|) supports the idea of the burst be- waf , ^ repoitgd fey Wekh ( mQ) afld confirmed afterwar ds by 

Send offprint requests to: A. Monreal-Ibero several works (e.g. |Walsh & Roy||1989|). 

* Based on observations collected at the European Organisation Narrow band imaging showed that as a consequence of its 

for Astronomical Research in the Southern Hemisphere, Chile (ESO starburst nature, the ionized gas in NGC 5253 pr esents a com- 

Programme 078.B-0043 and 383.B-0043). plex structure that includes filaments and arcs (e.g. Calzetti et al. 
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12004 . Therefore, proper characterization of the properties of 
the ionized gas in NGC 5253 would benefit from high qual- 
ity two-dimensional spectral mapping over a contiguous area 
of interest. Nowadays, Integral Field Spectroscopy (IFS) facil- 
ities represent the obvious choice to obtain this kind of infor- 
mation and several works devoted to the study of BCDs us- 
ing this approach have been published in recen t years (e.g. 
IzotovetalJ l2006at iGarcia-Lorenzo et al.l 120081: iJames etaU 
20091 IWestaoquette et al.ll2010t). ~_ 

In particular. iMonreal-Ibero eta! .] (1201 Obi hereafter Paper I) 
carried out a detailed study of the central region in NGC 525 3 
using IFS data collected with FLAMES dPasquiniet alj|2002l) . 
One of the topics studied in more detail there was the kinematics 
of the ionized gas. We found that the line profiles were complex, 
needing up to three components to reproduce them. In the main 
Giant Hn Region (GHuR, see Fig. Q]), one of them presented 
supersonic widths while the other two were relatively narrower. 
Moreover, the broad component presented a nd excess in nitrogen 
of ~1.4 times than that of the narrow one (iMonreal-ibero et all 
1201 lbl) . This was consistent with a scenario where the two SSCs 
produce an outflow that encounters the previously quiescent gas. 
Also, we delimited very precisely the area polluted with extra 
nitrogen, based on the excess of the [Nn]/i6584/Ho' ratio with 
respect to [S n]AA67 17,673 I /Ha and clearly demonstrated that, 
at least, some of the Wolf-Rayet (WR) star population cannot be 
the cause of this enrichment. Moreover, we could resolve a long- 
standing issue regardin g the elusive He n emission at A4686 A. 
Camp bell et al.l (U986T) mentioned a possible detection of He n 
but this result was not certainly confirmed afterwards. We de- 
tected several localized areas with clear He H/14686 detection al- 
though, rather puzzlingly, not coincident in general with the area 
exhibiting extra nitrogen. 

With the limited spectral range utilized in that work, we 
could only determine relative nitrogen abundances. Moreover, 
while interpreting the increase of [Nn]/l6584/[Sn]/L16717,6731 
as caused by an enhancement in nitrogen abundance is the most 
natural explan ation and is supported b y previous observational 
evidence (e.g. Kob ulnickv et alJ [l997). this is not the only al- 
ternative. Instead, specific combinations of ionization structure 
and metallicity gradient could also produce similar observables. 
Since dwarf galaxies are not expected to have strong metallicity 
gradients in a similar manner to large spirals, this option seems 
unlikely. However, it could not be rejected a priori in our previ- 
ous work. 

Her e, we present the natural continuation of the work in 
Paper I. To overcome the mentioned drawbacks, we obtained 
new FLAMES observations that allow us to map missing physi- 
cal properties such as electron temperature (T e ) and local degree 
of excitation. In this manner we are able to: i) map the chemi- 
cal content in the central part of the galaxy and detect inhomo- 
geneities (if any) in an unbiased manner; ii) evaluate how well 
relative abundance tracers based on strong emission lines repro- 
duce the values derived from direct measurements; iii) check if 
the derived physica l and ch emical structure is consistent with the 
picture sketched in PaperTl 

The paper is organized as follows: Sec. [2] describe the ob- 
servations and technical details about the data reduction and 
processing necessary to extract the required information for the 
analysis. Sec. [3] contains the results from one Gaussian fitting 
for the physical (i.e. T e , n e , ionization degree) and chemical (i.e. 
metallicity and relative abundances) properties both in 2D and 
in particularly interesting apertures. Sec. 2] w iH explore how the 
derived maps trace the 3D physical structure of the central part 
of NGC 5253, and the reliability (or limitations) of strong line 



based tracers to determine the relative abundance in nitrogen. 
Basic i nformation for NGC 5253 can be found in Table 1 in 

iPaper Ii 




Fig. 1. False colour image in filters F435W (B, blue chan- 
nel), F658N (Ha, green channel), and F814W (I, red chan- 
nel) for NGC 5253 using HST-ACS images (programme 10608, 
PLVacca). The area covered in ESO program 383.B-0043 
(078.B-0043) is marked with a white (grey) rectangle while the 
main GH nR of the galaxy is labeled in orange. The orientation 
and scale for a distance of 3.8 Mpc are indicated. 



2. Observations and data processing 

2.1. Observations 

Data were obtained with th e Fibre Large Array Multi Element 
Spectrograph, FLAMES (Pasu uini et alJ l2002h at Kueyen, 
Telescope Unit 2 of the 8 m VLT at ESO's observatory on 
Paranal. We used the ARGUS Integral Field Unit with the sam- 
pling of 0.52"/lens. This permits coverage of a field of view 
(fov) of 1 1'/5 x7'/3. Utilized gratings were L385.7 (LR1), L427.2 
(LR2), and L479.7 (LR3) and L682.2 (LR6). 

Data for the two last gratings were obtained in visitor mode 
on February 10, 2007 ( program me 078.B-0043). Details of these 
observations appear in IPaper fl In addition to these, data for the 
L385.7 (LR1) and L427.2 (LR2) were obtained in service mode 
during June 2009 (programme 383.B-0043). The spectral range, 
resolving power, exposure time and airmass for each configura- 
tion are listed in Table Q] Seeing ranged typically between 0'.'4 
and 172 with a median(+ standard deviation) of ~ 078 (+072). 
All the data were taken under clear conditions. In addition, stan- 
dard sets of calibration files were obtained. These included con- 
tinuum and ThAr arc lamps exposures as well as frames for the 
spectrophotometric standard stars HR 7596 and LTT 3218 for 
the L385.7 (LR1) and L427.2 (LR2) gratings respectively. 

Due to a guiding problem, there was a ~ 279 offset between 
the pointings on runs 078.B-0043 and 383.B-0043. In spite of 
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Table 1. Observation log 



Grating 


Spectral range 


Resolution 


texp 


Airmass 




(A) 




(s) 




L385.7 


3 610^1081 


12 800 


21 x 895 


1.00-1.10 


L427.2 


3 964-4 567 


10200 


9x895 


1.00-1.09 



it, the main emitting regions were covered by all the configura- 
tions. The precise area covered in each run is shown in Figure Q] 
which contains the FLAMES fov over-plotted on an HST B, Ha, 
I colour image. 

2.2. Data reduction 

The process ing of th e data on NGC 5253 for run 078.B-0043 is 
described in iPaper fl Those for run 383.B-0043 were processed 
using a combination of the pipeline provided by ESO (version 
2.8.913 y i a esorex, version 3.9.0 and some IRAF0 routines. 

The corresponding master calibration files were created with 
the ESO pipeline tasks gimasterbias, gimasterflat and 
giwavecalibration. In particular, fifteen bias frames were 
used for the masterbias creation while six and three continuum 
lamp exposures were needed for the creation of the L385.7 and 
L427.2 masterflats, respectively. In order to remove the cosmic 
rays in the target exposures, we combined all the individual 
frames for a given grating using imcombine with IRAF. After 
that, each combined frame was processed using the ESO pipeline 
in order to perform bias subtraction, spectral tracing and extrac- 
tion, wavelength calibration and correction of fibre transmission 
and data cube creation. 

In both settings, four isolated and evenly distributed arc lines 
were fitted by a Gaussian in each spectrum in order to estimate 
the uncertainties in the wavelength calibration as well as the in- 
strumental width. Centroids of the lines were determined with 
an accuracy of ~ 0.005 A which, for the covered spectral range, 
translates into velocities of ~0.3-0.4 km s _1 . We measured a 
spectral resolution of FWHM -0.329 A and 0.446 A for the 
L385.7 and L427.2, respectively. This translates to <Ti nstru - 10.9 
and 13.3 km s _1 . 

In the next step, the sky background was subtracted. For that 
purpose, we created a good signal-to-noise (S/N) spectrum by 
averaging the spectra of the sky fibres in each combined frame 
which was subtracted from every spectrum. Fibers suffering sig- 
nificant contamination due to cross-talk of adjacent object fibres 
were excluded from this combination. 

Relative flux calibration was performed within IRAF. The 
sky in the standard star frames was subtracted following the 
same methodology as with the science frames and then a spec- 
trum was formed for each calibration star by co-adding all the 
fibres of each standard star frame. Then, a sensitivity function 
was determined with the IRAF tasks standard and sensfunc 
and science frames were calibrated with calibrate. 

Finally, since data were taken in clear conditions, it was nec- 
essary to place both datacubes on the same (relative) flux scale. 
We scaled each spectrum in the LR1 mode to the same level as 



1 http://www.eso.org/projects/dfs/dfs-shared/web/vlt/vlt-instrament- 
pipelines.html. 

2 The Image Reduction and Analysis Facility IRAF is distributed by 
the National Optical Astronomy Observatories which is operated by the 
association of Universities for Research in Astronomy, Inc. under coop- 
erative agreement with the National Science Foundation. 



that of LR2 using the H7 and [Ne m]/l3967 emission lines, which 
were present in both instrumental set-ups. 

2.3. Subtracting the emission of the stellar population 

The spectral region analyzed in this work is rich in spectral fea- 
tures caused by the underlying stellar population (e.g. Balmer 
absorption lines). Their effect on the estimation of the flux in 
the emission lines for the gas is negligible for the brightest lines 
and/or those in the area of the bright GH nR. However, they may 
affect the measurements in the area of lower surface brightness 
and, especially in those spaxels associated with the relatively 
older clusters (-10 Myrfl associated to the peak of emission in 
the continuum jj3 (see IPaper J. for the definition of the peaks of 
emission). 

To estimate and correct this effect, we modelled and 
subtracted out the contribution of the underlying stellar 
continuum using the STARLIG HT^ spectral synthesis code 
(ICid Fernandes et ai1l2005l 12009). This code reproduces a given 
observed spectrum by selecting a linear combination of a sub- 
set of N* spectral components from a pre-defined set of base 
spectra. In our particular case, we utilized as base spect ra a set 
of single star populations from Br uzual & Ch ariot (20 0"3l). These 
are ba sed on the Padova 2000 evolutionary tracks dGirardi et alj 
2000) and assume Salpeter initial mass function between 0.1 
and 100 M H . Since NGC 52 53 has a metallicity of Z-0.3 
(e.g.[kobulnic kv et alJll999T) . we utilized only base spectra with 
Z =0.004 and 0.008. For each metallicity, we selected a set of 18 
spectra with ages ranging from 1 Myr to 2.5 Gyr. We allowed for 
a single extinction for all the base spectra that w as modeled as 
a unif orm dust screen with the extinction law by Cardelli et al. 
(fl989l 

The stellar spectral energy distribution at each spaxel was 
independently modelled for each cube. The spectral regions uti- 
lized in the fits were 3 650-4 060 A and the whole spectral range 
for the LR1 and LR2 gratings, respectively. All the main outputs 
for each spaxel were then reformatted in order to create three 
cubes per grating with the total, gaseous and modelled stellar 
emission. Additionally, the information associated to auxiliary 
properties such as stellar extinction in the V band, the stellar ve- 
locity and velocity dispersion, the reduced x 1 an d the absolute 
deviation of the fit (in %) were reformatted and saved as a file 
suitable for manipulation with standard astronomical software. 
Hereafter, we will use both terms map and image to refer to these 
kind of files. 

Representative examples of the achieved subtraction of the 
stellar component are shown in Fig. [2] which contains the ob- 
served, modeled and subtracted spectra for spaxels in the lowest 
surface brightness area, at the peak of emission in the contin- 
uum j)3 and in the peak of emission for the ionized gas. The 
subtraction of the stellar continuum is satisfactory in all three 
cases. Note that given the spectral resolution of the FLAMES 
data (larger than the one of the base spectra), some fine struc- 
ture (e.g. at -3 930 A) cannot be reproduced. However, this does 
not affect the recovery of the gaseous emission in which we are 
primarily interested. 



3 Ages for the clusters associated to this knot were incorrectly qu oted 
in our previous work. The ages estimated by iHarris et"aTl (120041) , for 
their clusters 3 and 5 were 8 and 1 1 Myr, respectively. 

4 http://starlight.ufsc.br/index.php?section= 1 

5 We assumed 12 + log(0/H) Q = 8.66 jAsplund et all2004h . 
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Fig. 2. Representative spectra to illustrate the modelling and subtraction of the underlying stellar population emission for the LR1 
(left) and LR2 (right) data cubes. The selected spaxels show typical emission in the lowest surface brightness area (first row), the 
peak of continuum emission jj3, associated to two relatively old clusters (second row), and the peak of emission in Hy (third row). 
Each individual graphic contains the observed, modelled and subtracted spectrum in black, cyan and orange, respectively. Below 
this main graphic there is an auxiliary one showing a zoom of the subtracted spectrum in black together with the spectral regions 
that were masked from the continuum fitting (green). 



2.4. Line fitting and map creation 

As a first step, information from the relevant emission lines for 
each individual spaxel was obtained by fitting Gaussian func- 
tions in a se mi-automatic way using the IDL based routine 
MPFITEXPR (Ma rkwardlll2009l) . Most of the lines (i.e. Hy- 
H8, [Om]/l4363, HeiA4009, He 1^4026, He l!4388, He l!4471, 
[Sii]/14074 and [Sn]/14076) were strong enough over most 
of the fov to be fitted independently. Exceptional cases were 
[Ne iii]/13967 and He 1/13965 that were fitted simultaneously with 
H7, to assure a proper deblending of the emission lines and 
[On]/l3727 and [On]/i3729 that were fitted simultaneously. For 
the rest of lines, the kinematic results forH7, present in both LR1 
and LR2, were utilized as additional constrains. Several (strong) 



lines showed signs of asymmetries and/or multiple components 
in their profiles for a large number of spaxels. When possible, 
multi-component fits to these lines were also performed trying 
to assure a continuity in the derived physical (and chemical) 
properties between adjacent spaxels. In the present work, this 
multi-component fit will only be utilized in the discussion of the 
electron density structure. Finally, as we did with the fitting of 
the underlying stellar population, we used the derived quantities 
together with the position within the data-cube for each spaxel 
to create a map. 

On account of the different pointings between the 078. B- 
0043 and 382.B-0043 observing runs, an important issue for this 
analysis was assuring a proper match between maps for quan- 
tities derived from the LR3 and LR6 cubes and those of the 
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Ax (arcsec) 



Fig. 3. Ionized gas distribution as traced by the [O n] emission 
line doublet. We have over-plotted contours corresponding to the 
HST-ACS images in the F435W filter (programme 10609, P.I.: 
Vacca) convolved with a Gaussian of 0V8 to simulate the seeing 
at Paranal. The position of the three main peaks of continuum 
emission are marked with crosses. The map is presented in log- 
arithmic scale in order to emphasize the relevant morphological 
features and cover a range of 1 .4 dex. Flux units are arbitrary. 
Utilized apertures to extract the selected spectra have been out- 
lined according to the following color code: Cyan - Knot fll; 
Blue - Knot j}2; Red - Knot f)3. Note the existence of three dead 
spaxels at ~ [3'.'. 5, l'.'O] as well as absence of signal in the spax- 
els at the four comers of the field of view. All of these spaxels 
will be marked hereafter by white rectangles. The scale and ori- 
entation are indicated. 



LR1 and LR2. As a first guess we used the three main peaks 
in the continuum images as our reference to align the images. 
However, given the relatively small fov of the IFU, the extended 
nature of knots |jl and ((3 and that knot (J2 falls at the edge of 
the IFU in the LR1 and LR2 data, the derived offset were not 
accurate enough. To refine this first guess, we created maps for 
ratios involving lines observed with the two different pointings 
(e.g. [Oiii]/15007/[Oiii]4363) using a grid of offsets centered in 
our first guess and in steps of 0.5 spaxels (=0'/26~l/4th seeing). 
We consider as final offset [2.5 spa,5.0 spa] (=[1'.'3,2'.'6]) which 
was the one providing the best continuity in adjacent spaxels for 
the variation of the line ratios. If the map of a given quantity 
was based on information from both runs, when possible, we 
independently carried out the calculations for 078.B-0043 and 
383.B-0043. Only at the end, the offset was applied to the infor- 
mation associated to the 078.B-0043 run before carrying out the 
final calculation involving data from both runs. In this manner, 
the number of spatial interpolations was minimized. 

3. Results 

3.1. Gaseous emission in selected apertures 

The main aim of the present work is to provide a full 2D char- 
acterization of the physical and chemical properties of the ion- 
ized gas in the central region of the galaxy. However, there are 
reasons to perform a more detailed analysis in specific regions. 
Firstly, we can define very precisely the aperture utilized to ex- 
tract a given spectrum, avoiding the limitations associated to the 
positioning of a slit and assuring that the spectrum is spatially 
associated with a given feature of interest (e.g. the location of a 



Table 2. Observed fluxes for lines in the LR1 and the LR2 grat- 
ings with respect to F(H/3). (See text for details of the scaling of 
the spectra ). 



Jo 


Ion 


fW 




F(i)/F(H/3) 




(A) 






Knot ((I 


Knot j}2 


Knot (|3 


3669.47 


Hi 


0.268 


0.17±0.05 


0.42±0.05 




3671.48 


Hi 


0.268 


0.24±0.05 


0.24±0.06 




3673.76 


Hi 


0.268 


0.35±0.05 


0.46±0.06 




3676.37 


Hi 


0.265 


0.53±0.06 


0.57±0.07 




3679.36 


Hi 


0.265 


0.52±0.04 


0.69±0.11 




3682.81 


Hi 


0.265 


0.67±0.05 


0.77±0.08 




3686.83 


Hi 


0.263 


0.75±0.05 


0.99±0.08 




3691.56 


Hi 


0.263 


0.83±0.06 


0.91±0.06 




3697.15 


Hi 


0.261 


1.10±0.08 


0.91±0.06 




3703.86 


Hi 


0.261 


1.30±0.08 


1.25±0.06 




3705.04 


He i 


0.261 


0.70±0.07 


0.57±0.06 




3711.97 


Hi 


0.259 


1.37±0.08 


1.71±0.09 




3721.83 


[Shi] 


0.257 


2.95±0.16 


2.68±0.09 




3726.03 


[On] 


0.255 


43.82±2.73 


75.36±2.05 


101.88±3.97 


3728.82 


[On] 


0.255 


46.35±2.77 


93.87±2.61 


144.38±6.84 


3734. 17 


H i 


ft OK 1 ^ 


1 GG-i-ft 1 A 


1 QG_i_n ftQ 

i .vy±u.us 


1 .zo±u.zu 


3750.15 


Hi 


0.251 


2.70±0.16 


2.80±0.10 


3.24±0.40 


3770.63 


Hi 


0.247 


3.34±0.18 


3.54±0.09 


2.51±0.32 


3797.90 


Hi 


0.241 


4.59±0.25 


4.69±0.12 


4.56±0.33 


3819.61 


He i 


0.237 


0.86±0.06 


0.88±0.04 


0.65±0.20 


3835.39 


Hi 


0.235 


6.61±0.35 


6.74±0.15 


6.65±0.48 


3868.75 


[Nem] 


0.227 


48.87±3.10 


31.73±0.93 


22.34±0.89 


3889.05 


H i+He i 


0.223 


16.52±0.90 


18.1 1±0.43 


15.91±0.67 


3964.73 


He i 


0.209 


0.57±0.35 


0.59±0.19 


0.77±0.45 


3967.46 


[Ne m] 


0.207 


15.88±0.82 


10.38±0.31 


8.77±0.59 


3970.07 


Hi 


0.207 


15.80±0.82 


15.80±0.38 


15.80±0.72 


4009.22 


He i 


0.198 


0.20±0.03<*> 






4026.21 


He i 


0.194 


1.73±0.09 


1.69±0.05 


0.81±0.19 


4068.60 


[S ii] 


0.187 


1.40±0.08 


1.41 ±0.05 


3.28±0.24 


4076.35 


[Sn] 


0.184 


0.46±0.03 


0.52±0.03 


1.25±0.16 


4097.26 


On 


0.180 


0.05±0.01 






4101.74 


Hi 


0.180 


25.79±1.37 


25.42±0.58 


25.19±0.89 


4120.82 


Hei 


0.176 


0.20±0.02 


0.09±0.03 




4143.76 


He i 


0.172 


0.23±0.02 


0.16±0.02 




4243.97 


[Feu] 


0.153 


0.03±0.01 


0.05±0.01 




4267.15 


Cn 


0.147 


0.06±0.01 


0.10±0.01 




4287.40 


[Feu] 


0.144 


0.12±0.01 


0.24±0.02 


1.23±0.15 


4340.47 


Hi 


0.133 


46.60±2.87 


46.60±1.12 


46.60±1.71 


4359.34 


[Feu] 


0.129 


0.06±0.01 


0.13±0.01 


0.53±0.10 


4363.21 


[Om] 


0.129 


7.77±0.39 


3.70±0.08 


2.64±0.21 


4368.25 


Oi 


0.127 


0.03±0.01 






4387.93 


He i 


0.122 


0.49±0.03 


0.45±0.02 




4413.78 


[Fen] 


0.118 


0.07±0.01 


0.01±0.01 




4416.27 


[Feu] 


0.115 


0.04±0.01 


0.03±0.01 




4437.55 


He i 


0.110 


4.48±0.25 






4452.11 


[Feu] 


0.108 


0.13±0.01 






4471.48 


He i 


0.103 


0.08±0.01 


0.19±0.01 


0.36±0.18 



This line presented a "secondary peak" towards the blue which cor- 
respond to a weak [Fern] line at /14008. At the resolution of the 
present data a proper deblending of the lines was not feasible. 



star cluster). As an example, Lop ez-Sanchez et al 1 (12007b char- 
acterized the properties of our knot (j2 (their UV-1). However, 
neither their positions HII-1 nor HII-2 sample the gas associ- 
ated to the peak of emission in Ha, but some emission asso- 
ciated to the GHuR towards the west and south of the peak 
of emission. Secondly, a comparison of the measurements pro- 
vided here, with those already existing in the literature serves as 
a sanity check which will reinforce the results obtained in the 
2D mapping. Thirdly, most of the present methodologies used 
in the analysis of the physical and chemical properties of the 
ionized gas were developed to study complete H n regions. At 
present, the astronomical community is still in the process of un- 
derstanding under wh ich conditions these can be applied to small 
portio ns of them (e.g. Perez -Mon tero et aT1l201 lUErcolano et al.l 
2012). The areas analyzed in this section will be closer to this 
notion of complete H n than individual spaxels. Finally, by com- 
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Table 3. Observed fluxes for lines in the LR3 and the LR6 grat- 
ings with respect to F(H/3). (See text for details of the scaling of 
the spectra ). 



M 


Ion 


ft y\ 
JW 




t{A)/t{tip) 




(A) 






Knot pi 


Knot p, 


Knot p3 


4562.60 


Mg i] 


0.080 


0. 1 1 ±0.01 


0.19±0.01 


0.86± 0.19 


4571.00 


Mg i] 


0.078 


0.12±0.01 


0.16±0.01 


0.74± 0.17 


4649.13 


On 


0.057 


0.08±0.01 






4658.10 


[Fe m] 


0.054 


0,96±0.05 


0.80±0.03 


2.34± 0.19 


4701.53 


[Fe hi] 


0.043 


0.27 ±0.01 


0.23±0.01 


0.84± 0.19 


4711.37 


[Ariv] 


0.041 


1.17±0.06 


0.19±0.00 


0.57± 0.02 


4713.14 


He i 


0.041 


0.65±0.05 


0.41±0.01 


0.19± 0.16 


4740.16 


[ Ar iv] 


0.032 


1 . 1 8±0.05 


0.17±0.01 




4754.83 


[ Fe 1 1 1 1 


0.030 


u. iy±u.ui 


A IO i fl 111 




4769.60 


[Fe m] 


0.024 


0. 10±0.01 


0.10±0.01 




4814.55 


[Fe n] 


0.013 


0.06±0.01 






4861.33 


H i 


0.000 


100.00±4.36 


100.00±1.82 


1 f\t\ t\f\ i 1 CI 

100. 00± 2.51 


4881.00 


[Fe hi] 


-0.005 


0.37±0.02 


0.25 ±0.02 




hah r\i 

4921.93 


He i 


—0.016 


1 1 a i <~\ r\c 
1.14±0.05 


1.03±0.03 




4931.32 


[O m] 


-0.019 


0. 16±0.01 


O.lo±0.02 




4958.91 


[O m] 


-0.027 


232.56±10.88 


162.83±3.25 


1 14.29± 3.25 


4985.90 


[Fe ni] 


-0.032 


0.47±0.03 


0.64±0.03 


2.78±0.26 


5006.84 


[O in] 


-0.040 


732.56±34.05 


508.85±10.23 


358. 10± 10.25 


5015.68 


He i 


-0.043 


1.87±0.09 


2.19±0.06 


2. 11 ±0.25 


5041.03 


Si ii 


-0.049 


0.24±0.02 


0.16±0.01 




5055.98 


Si ii 


-0.055 




0.16±0.01 




6548.03 


[Nil] 


-0.312 


10.7()±0.46 


5.88±0.16 


11.33± 0.29 


6562.82 


Hi 


-0.314 


400.00± 16.86 


391.15±7.37 


415.24± 10.18 


6583.41 


[Nil] 


-0.316 


32.67±1.39 


17.96±0.49 


35.05± 0.80 


6678.15 


He i 


-0.329 


4.55±0.20 


4.20±0.08 




6716.47 


[S ii] 


-0.332 


13.95±0.60 


20.53±0.41 


55.24±1.32 


6730.85 


[Sn] 


-0.335 


12.91±0.55 


16.55±0.34 


40.76±0.98 


7002.23 


Oi 


-0.365 


0.09±0.01 


0.10±0.01 




7065.28 


He i 


-0.375 


7.73±0.30 


3.85±0.08 


2.83±0.19 


7135.78 


[Ar m] 


-0.384 


16.34±0.73 


12.65±0.26 


10.67±0.30 



Table 4. Physical properties of the ionized gas. 





Knotftl 


Knot p. 


Knot (t3 


E(B - V) 


0.31±0.04 


0.29+0.02 


0.35+0.04 


^(cnr^aOn]) 
n e (cnr 3 )([Sn]) 
«,(cm" 3 ) ([Ariv]) 
n c (cnr 3 ) ([Fern]) 


310 +15 ° 

-120 

4 600+ 1200 

^ -30 


140 +5 ° 
170+^ 
2 8001™ 
155+20 


ap+100 

-60 

70+30 

'"-25 
7 -20 


r e (K)([Om]) 
r c (K)([On]) 
r e (K)([Sn]) 


11 570+go 

11360 

97 Oo+i70o 
v ivv _ 13(x) 


10280 + | 3 ° 

10470 
8 700 + ^» 


10 320+ 350 

—300 

10 500 

8 900+ 870 



- Ratios for n e : [On]: [O ii],13727/[0 n]/13729; [Sn]: 
[Sii]/16717/[Sn]/16731; [Ariv]: [Ariv],14711/[Ariv]/14740; [Fern]: 
[Fem]/14986/[Fem]/l4658 

- Ratios for T e : [Om]: [Om],U4959,5007/[O m]/14363; [Sn]: 

[S n]^6716,6731/[S n]^4069,4076. Note that r c ([On]) was not 
obtained directly from measurements of [O n] lines. Instead, it was 
derived fro m r,([Om]) an d n e ([O n]) according to the expressions 
provided in lGarcia- Benito (2009). 



bining the signal of a collection of spaxels, a larger signal-to- 
noise (S/N) ratio is achieved and fainter lines can be detected. 

As we saw in Paper I, the stellar emission is dominated by 
three peaks of emission named as knots jj2 and (J3 (see Tab. 3 
of Paper I to establish the correspondence between this nomen- 
clature and those of previous studies). Note that knot #1 is associ- 
ated to the main GH nR marked in Fig.Q] The spectra associated 
to them have been co-added and extracted in order to determine 
relative line intensities in their surrounding ionized gas . 

Utilized apertures in the LR1 and LR2 data are outlined in 
Fig. [3] For the LR3 and LR6, we applied an offset of 2 and 5 
spaxels in the x and y direction respectively. Since the offset be- 



Table 5. Ionic chemical abundances of the ionized gas. 





Knotftl 


Knot tt2 


Knot (t3 


12 + log(0 /H ) 


7.33±0.10 


7.73±0.04 


"7 O O . f\ f\C\ 

7.88+0.09 


12 + log(0 /rt ) 


8.18±0.09 


o 1 "7 . n r\i 

8.17±0.03 


8.03+0.07 


12 + log((9///) 


8.24±0.09 


8.32±0.04 


8.26+0.08 


12 + log(A /H ) 


5.7l±0.07 


5.96±0.03 


6.31+0.05 


log(S + /0 + ) 


-l.62±0.l7 


-1.77±0.07 


-1.57+0.14 


12 + log(/V + /# + ) 


6.52±0.07 


6.36±0.03 


6.62+0.06 


12 + log(/V///) 


7.43±0.07 


6.94±0.03 


7.00+0.06 


lL.t (Ar ) 


8.04 


3.84 


2.40 


lOg(N/U) 


— U.Sl+U.Zj 


i T7_i_n 1 1 
— l.j l +U. 1 1 


— 1 .ZO+U.Zz 


12 + log(/Ve 2+ /// + ) 


7.49±0.10 


7.50±0.04 


7.34+0.09 


ICF (Ne) 


1.14 


1.35 


1.71 


12 + log(Ne/H) 


7.55±0.10 


7.63±0.04 


7.57±0.09 


log(Ne/0) 


-0.69±0.28 


-0.69±0.12 


-0.69±0.25 


12 + log(Ar 2+ /# + ) 


5.89±0.07 


5.90±0.03 


5.80+0.05 


12 + log(Ar 3+ /H + ) 


4.84±0.08 


4.15±0.05 




ICF (Ar) w 


1.02±0.04 


1.07±0.05 


1.45+1.91 


12 + log(Ar/H) 


5.90±0.22 


5.93±0.12 


5.96±1.45 


log(Ar/0) 


-2.34±0.30 


-2.38±0.16 


-2.30+1.53 


12 + log(Fe 2+ /# + ) 


5.58±0.09 


5.66±0.04 


6.12+0.09 



w ICF (Ar 2+ )for knot (t 3. 



tween the two pointings was determined as [2.5 spa,5.0 spa], this 
implies a small difference in the selected apertures from the old 
and the new data which introduces an extra uncertainty when 
comparing fluxes of lines belonging to different set of data. In 
consequence, the line fluxes for the LR3 and the LR6 gratings 
were directly measured with respect to H/J and Ha while those 
for the LR1 and LR2 gratings where measured with respect to 
Hy and H7, and then w e assumed the theorietical Balmer line 
intensities obtained from lStorev & Hummer! (1 1995b for Case B, 
T e = 10 4 K and n e = 100 cm 4 . Most of the lines were mea- 
sured independently by fitting the emission to a single Gaussian 
profile using rapfit. Exceptionally, lines in the pairs H16 
and He^3705, [Fen]^4414 and [Fen]^4416, and [Ariv]/I4711 
and He 1/147 13 were simultaneously fitted due to their prox- 
imity in wavelength. Similarly, all the three lines He 1/13965, 
[Nem]/13967 and H7 were fitted at once. Measurements for all 
detected lines are compiled in Tables |2]and[3] 

3.1 .1 . Physical conditions of the ionized gas 

Derived physical conditions for the ionized gas are listed in 
Table [4] The values of reddening were derived from the mea- 
sured Hor/H/J line ratio following the methodology described in 
Paper I. Electron temperature (T e ) and density (n e ) based on 
[Sn], [On] and [Om] emission lines (and thus tracing differ- 
ent layers in the ionizati on structure) were d erived using the ex- 
pressions provided by Garcf a-Benitol ([20090- For that we pro- 
ceeded as follows: we assumed an initial n e of 100 cirT 3 to ob- 
tain a first guess of the different T e 's. The resulting r c ([On]) 
and r c ([Sn]) were used as input to obtain new estimates of 
« e ([OffJ) and n e ([Sn]) and the process was iterated until con- 

6 http://www.iac.es/consolider-ingenio- 
gtc/index.php?option=com_content&view=article&id=223:qa- 
spatially-resolved-study-of-ionized-regions-in-galaxies-at-different- 
scales&catid=45 :tesis&Itemid=65 



6 



A. Monreal-Ibero et al.: Physical and chemical conditions of the ionized gas in NGC 5253 



vergence. Typically one or two iterations were enough. These 
functions reproduce the predictions of the t ask temden, based on 
the f ivel program (Sha w & Dufourll995l) included in the IRAF 
p ackage nebular. They used th e same atomic coefficients as 
in lPerez-Montero & Diazl d2003l) . exce pt for Q + for w hich they 
used the transition prob abilities f rom | Zeippenl (1 19821) and the 
collision strengths from iPradhanl (119761) . r.dOnll was not in- 
dependently derived from our data, since the [O ii]/1/17320,7330 
doublet was not covered in these observations. Instead, we de- 
rived r g ([O nl) from T,,(\Om]) and »,,([() nl) using the models 
presented bv lPerez-Montero & Diazl(l2003l) . A s an improvemen t 
to former relations based on modelling (e.g. IStasinskal fl990). 
these take into account depende ncies on n P . At low densities , 
for a given r<,([Om]), models bv lPerez-Montero"& Dfazl(l200l 
predict on l y sligh tly (i.e. <2%) lower temperatures than those by 
IStasinskal d 19901) while in the densest zones of NGC 5253 (e.g. 
our knot jjl), estimated r f ([On]) can be up to ~ 17% smaller. 
For « e (fFe in]), we used the line ratios tabulated bv lKeenan et alj 
(1200 ll) . assuming T e = 12 000 K. Note that for the utilized emis- 
sion lines, the dependence of the derived n e on the assumed tem- 
perature is negligible. Finally, n e ,([Ariv]) was derived using di- 
rectly temden and assuming r e ([Ariv])=7' f ,([Oiii]). 

Derived T e and n e for our knots jjl and fi2 agree well w ithin 
the uncertainties with the values reported by ISidolil (1201 Ol) for 
his apertures A and B, respectively. Addit ional measurements 
for T e and n e in knot |j2 are pr ovided by Lopez-Sanchez etaU 
d2007l) . lKobulnickv et alJ 1 19971) and lGuseva et alJ d201 ll) . While 
our derived T e 's are in relatively good agreement, with T e ([0 ra]) 
and r<,([Sn]) being slightly lower and larger respectively than 
the values reported in these works, there is a discrepancy be- 
tween the different measurements of n e in the literature. A com- 
parison of ou r n e (\S n]) for knot jj 1 and the values reported by 
iGuseva et alJ (1201 ll) for their apertures CI and P2 shows a dis- 
crepancy between the reported values. Since n e varies a lot in 
this area (a factor of ~2-3 on scales of ~(X'5, see Fig.[9]in Sec. 
14.1b . the precise definition of the aperture (i.e. size, position and 
shape) in each case seems the most plausible explanation for 
this discrepancy. Our n e ~ 150 cm in knot jj2 is consistent 
within the errors with being below the low density limit as re- 
ported by iKobulnickv et al.l d 19971). However, it is a factor ~ 2 
lower than the values reported bv ILopez-Sanchez et al.l (120071) . 
While estimated uncertainties could account for the differences 
in n e ([S n]), one must resort to other causes, such as differences 
in the definition of the aperture, to explain the discrepancies in 
«e([On]). 

3.1.2. Chemical abundances 

Abundances for the different ionic species are listed in Table 
an d were derived using the relations provided in Appendix B.2 
of iGarcfa-Benito (2009). These are appropriate fittin gs to the 
results of the IRAF task ionic (Shaw & Dufo urj|l995l). based 
on the 5-level atom program developed by De Robertis et alJ 
1987), and follow the functional form given by iPagel et all 
1992). For the abundances of neon, argon, iron and ++ , we 



utilized r f ([Om]). For + and nitrogen, we utilized r c ([On]), 
derived from r<,([Om]) and « e ([On]). Finally, for sulfur, we 
took r e ([Sn]). Differences in abundances due to the assumed 
T e are discussed below. To derive the total abundances, unseen 
ionization stages of a given element were taken into account 
by including the appropriate ionization factor (ICF) for each 
species when necessary, follo wing the prescriptions provided by 
Kingsburgh & Barlow ( 1994) for all the ele ments but argon. Fo r 
this element, we utilized those provided bv llzotov et al.l (Il994l) . 



Oxygen: We assumed O/H = (<9 + + 2+ )/H + . The derived 
values are in good agreement with those provided in the liter- 
ature using long- slit data (e.g. IGuseva et al.ll201 ll [Sidolill2010fc 
ILopez-Sanchez et al.ll2007l IKobulnickv et alJll99% . The largest 
differences are for knot fjl t hat shows abundan ces ~0.1 dex 
higher than those measured bv lGuseva et al. l d20TTl) . This is how- 
ever within the uncertainties. 

Sulfur: Deriv ed values are typically ~0.2 dex larger than 
those provided by Lope z-Sanchez et ail (|2007|). Differences can 
be attributed to the adopted ^.Lopez-Sanc hez et al.l d2007l) use 
the average of r e ([Sn]), r e ([On]), and r e ([Nn]) that can be 
~2 000 K larger than T e ([S n]). If a T e ([On]), derived from the 
n e ([On]) and r<,([Oin]), were utilized, our abundances would 
become typically ~0.15 dex smaller in agreement with the val- 
ues reported by these authors. A similar effect i s obser ved when 
comparing with the results reported by ISidolil (12010) who uti- 
liz ed TMOm]) to deriv e all the abundances and those reported 
bv IGuseva etai] d20lTl) who used a r e ([On]) -1500-2000 K 
larger than the T e {[S n]) utilized here. 

Nitrogen: We assumed N/O = N + /0 + , which is an accu- 
rate approximation to about ±20 % for nebulae with metallici- 
ties smaller than that in the LMC (lGarnettl|l990T) . In agreement 
with previous work, while knots (j2 and jj3 present relative abun- 
dances in nitrogen within the range ex pected for galaxies at this 
metallicity (see e.g. Mo lla et al.l2006l) . the aperture associated to 
knot jjl present a clear excess in nitrogen. Interestingly, although 
consistent with the N/ O expected for galaxies at this metallicity, 
knot jj3 has a slightly larger abundance than knot jj2. 

Neon: Since the ionization structure is similar to the one 
of the oxygen, we assumed Ne 2+ /Ne = 2+ /O. Reported val- 
ues agree wit hin the errors with previous measurements in 
the literature (IGuseva et al.l 120111: ILopez-Sanchez et al" 2007; 
IKobulnickv etalJI 19971) . 

Argon: Total argon abundance was de rived using t h e ICF s 
based on the expressions provided by llzotov etall dj994). 
Derive d ratios are within t he range of those previous l y pub- 
lished dGusevaetalJl201lt ILopez-Sanchez et all 120071: ISidolil 
1201 Ol) and consistent with a homogeneous Ar/O ratio across the 
face of the galaxy. 

Iron: Reported values f or Fe 2 + agree with those given by 



ILopez-Sanchez et all d2007l) 

To recap, abundances of the different ionic species in 
knot jjl and knot jj2 agree in general with those previously 
reporte d (Kobulnickv et al. 1997; ILopez-Sanchez et aD [2007; 
ISidolill2010ir Abundances in knot jj3 (not reported so far) are 
similar to those in knot jjl and knot jj2. With the exception of 
the relative abundance for N/O in knot jjl and maybe knot f)3, 
no chemical species seems overabundant in any of the selected 
apertures. This is in agreement with the relative enrichment in 
nitrogen previously widely reported. 

3.2. Mapping electron density and temperature tracers 

In iPaper H we found an « e ([Sn]) gradient declining from the 
peak of emission in Ha outwards and traced the gas at the highest 
densities by means of the n e ([Ar iv]). Also, our multi-component 
analysis of the [S n] line profiles tentatively suggested similar 
densities over the whole face of the main H n region for the broad 
component, while the narrow component presented somewhat 
lower (higher) ratios towards the NW (SE) part of the region. 
Given the restricted spectral range utilized in that work, the T e 
structure could not be derived. 

Here, we present maps for two additional tracers of n e : 
[On]/l3726/[On]^3729 and [Fera]^4986/[Fera]/14658. In com- 
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[OII]X3726/[OII]?l3729 




6 R 



[01 1 l]U4959,5007/[OII l]X4343 



2 0-2-4 -6 

Ax (arcsec) 

[Fel I l]A.4986/[Fel I l]X4658 
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CO 

o 

i— 

< 







- 1 1 1 1 1 1 




i i i 





-2 -4 -6 

Ax (arcsec) 



+1.80 

+1.42 
+1.05 
+0.68 
+0.30 



Fig. 4. Maps for the new line ratio sensitive to the electron den- 
sity. The position of the three main peaks of continuum emission 
are shown as crosses for reference. 



parison to [Sn]/l6717/[Sn]/l6731, the first ratio is sensitive 
to a slightly larger range o f densities, tracing to lower n e 
dOsterbroc k & Ferland 2006). Moreover, the [O n] lines are 
stronger than those of [S n]. Therefore, we will be able to more 
systematically discuss the differences in n e structure between 
different kinematic components, as suggested by the [S n] line 
ratio. The second ratio is sensitive to an e ven larger range of 
densities, from ~ 10 2 cm~ 3 to <10 7 cnr 3 (Kee nan et al] l2001) 
and has already been used with success in s imilar works to thi s 
one to reveal very high density locations dJames et al . 2009). 
Also, we will introduce for the first time maps for two trac- 
ers of T e in this galaxy: [Om]/U4959,5007/[Om]/l4343 and 
[Sn]/U6717,6731//U4069,4076. By comparing the maps in n e 
and T e according to the different tracers, we will be able to see 
how the physical properties of the gas are structured in the dif- 
ferent ionization layers (see Sec.|4]i. 

The maps for our new n e tracers are presented in Fig. |4] For 
that involving the [O n] emission lines, large values, associated 
to relatively large n e , are found in the main GH nR. There is a 
secondary maximum at the location of knot jj2 and ratios typical 
of low densities elsewhere. The map for the [Fe in] line ratios 
presents the inverse tendencies: the maximun of [O n] line ratio 
corresponds to the minimum of the [Fe m] line ratio. The two 
overal l structur es are similar to the one depicted by the [S n] line 
ratios dPaper J) . 

More interestingly, Fig. [5]presents the maps for the line ra- 
tios tracing the T e . Both of them display similar structure, with 
the absolute minimum at the location of the two SSCs (i.e. knot 
|}1), corresponding to the highest T e , ratios remaining relatively 
low in the GH uRand then tending to higher values in the rest of 
the fov. As with the selected apertures in Sec. I3.1.1I on each in- 
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Fig. 5. Maps for line ratio sensitive to T e tracing the high {up) 
and the low (bottom) ionization regime. The position of the three 
main peaks of continuum emission are shown as crosses for ref- 
erence. 



dividual spaxel the predicted T e from the [O in] line ratio differs 
from the one using the [S n] lines. We will explore in Sec.|4]how 
this is related to the ionization structure of the galaxy. 



3.3. Mapping tracers of the local ionization degree 

The relation between the [N n]/l6584/Ha and 
[S ii]/L167 17,673 1/Hor line ratios in galaxies is not mono- 
tonic, even though the S + and ionization potentials and 
critical densities are in a similar range. Instead, different 
ionization mechanisms (shocks, star formation, AGNs, e.g. 
Mon real-Ibero et al . 2006, 2010a) and/or physical and chemical 
conditions of the gas (metallicity, relative abundances and de- 
gree of ionization) trace different loci in th e [Nn]/16584/Ha vs. 
[S n]/L167 17,673 1/Ha diagram. In [Paper J, we attributed the dif- 
ferent structure in the [N n]/l6584/Ha and [S n]/U6717,6731/Ha 
maps to an excess in nitrogen. This was supported by evi- 
dence of star-formation as the main ionization mechanism and 
previous estimations of the ionization para meter at specific 
locati ons based on long-slit measurements ( Kobulnickv et al.l 
1997). However, certain combinations of specific ionization 
structures and metallicity gradients could reproduce a similar 
locus in the [Nn]i6584/Ha vs. [S n],U6717,6731/Ha diagram. 
This possibility can be ruled out with the present data. Here, 
we present the ionization structure in NGC 5253 as traced by 
different line ratios, while the next section will contain maps for 
the abundances of different heavy elements (oxygen, nitrogen, 
neon and argon). 
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Fig. 6. Maps for two tracers of the ionization degree. Top: 
[Om]/[On]. Bottom: [On]/H7. The position of the three main 
peaks of continuum emission are shown as crosses for reference. 



The degree of ionization can be traced by means of ratios 
of lines of the same element tracing two different ionization 
states. The upper panel of Fig. [6] contains a map for one of 
these ratios: [Om]/U4959,5007/[On]/U3726,3729. Note that in 
order to minimize uncertainties associated to aperture match- 
ing, absolute flux calibration and extinction, the [O in] and [O n] 
lines were measured with respect to HB and H7. Then, we as- 
sumend the theoretical B almer line intensities obtained from 
IStorev& Hummed (1 1995b for Case B, T e = 10 4 K and n e = 100 
cm 3 . The map shows that the ionization structure reproduces the 
morphology observed for the ionized gas. That is: i) line ratios 
tracing the highest ionization degree are associated to knot |Jl; 
ii) relatively high ionization is found in the main GH nR; iii) a 
secondary peak of high ionization is found around knot j}2; iv) 
low ionization degree is found in the rest of the fov, where the 
diffuse component of the ionized gas becomes more relevant. 

Also, if the metallicity is known, line ratios like 
[On]^3726,3729/H7 or [S n]/U67 17,673 1/Hor can also be used 
to trace the ionization degree. The map [On]/l/l3726,3729/H7 
is presented in the lower panel of Fig. [6] while that of 
[S u]AA67 17,673 l/Ha was included in PaperJ The similar 
structure in all the three maps suggests a lack of metallicity gra- 
dient in the galaxy as will be shown in the next section. 

In order to explore whether these three tracers predict con- 
sistent ionization degree for a given position, we show in 
Fig. [7] the relation between the different ratios for each in- 
dividual spaxel, using [Om]/U4959,5007/[On]^3726,3729 as 
reference. Also, we overplotted the locus of equal estimated 
ionization parame ter according to the relations proposed by 
iDfaz et al .] (120001) and assuming Z = 0.3 Z with a black 
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Fig. 7. Diagrams involving two ratios sensitive to the 
ionization degree. Top: [Om]X14959,5007/[Oii],U3726,3729 
vs. [On]/L13726,3729/H J e. Bottom: [Sn]/U6717,6731/Ha vs. 
[Om]X14959,5007/[Oii],M3726,3729). The locus of equal esti- 
m ated ionizat i on par ameter according to the relations proposed 
bv lDiaz et alj d2000t) and assuming Z = 0.3 Z is indicated with 
a black line. Only those spaxels with estimated uncertainties 
for the involved ratios smaller than 0.5 dex have been taken 
into account. Spaxels were grouped in three bins according their 
n e ([0 n]) as indicated in the upper right-hand corner of the indi- 
vidual diagrams. The total number of considered data points is 
indicated in the lower right-hand corner. 



line. For a given location, [O m]/U4959,5007/[On] ^3726,3729 
and [On]/l/l3726,2729/HjS predict a similar degree of ioniza- 
tion while the predictions of [S n)AA61 17,673 l/Ha would cor- 
respond to smaller ionization parameters. A similar result was 
found in a detailed analysi s of NGC 588, a GHnR in M 33 
dMonreal-Ibero et all 1201 lah . Both, the spatial variations of the 
different line ratios (Fig. |6]l and the observed excess in the 
[S n]/U67 17,673 l/Ha ratio when compared with photoioniza- 
tion models (Fig. [7} can be jointly explained as a 3D view of 
the ionization structure of the galaxy. Specifically, for a given 
spaxel (i.e. a given line of sight), the lower ionization species 
(e.g. S + ), delineate the more extended diffuse component while 
++ will be confined to different high ionization zones closer to 
the ionizing sources. Alternatively, the difference between the ra- 
tios found here and those predicted by the Diaz et al. models can 
be attributed to differences betwe en the relative S /O abundances 
in NG C 5253 (log(5/0) ~ - 1 .47. lKobulnickv et alJl997l: rSidoli 
2010) and those assumed in the models (log{S/0) = -1.71, 
iGrevesse & An ders 1989). A detailed modelling of the ioniza- 
tion structure of the galaxy - out of the scope of this work - could 
help to disentangle these two possibilities. 
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Fig. 8. Metallicity and relative abundances maps. From top to bottom and left to right: Oxygen abundance, and nitrogen, neon and 
argon relative abundances with respect to oxygen. In order to make easier the comparison between the different maps, all of them 
cover a range of 1 dex. 



3.4. Mapping abundances of heavy elements 



Following the same methodology as in Sec. 13.1.21 we derived 
maps for the abundances of oxygen, nitrogen, neon and argon. 
These are presented in Fig. [8] and show, as expected, a ho- 
mogeneous distribution within the uncertainties in all the el- 
ements but nitrogen. Note that the spaxel-to-spaxel variations 
in these maps are not dominated by the quality of the data 
(e.g. S/N ratio), but are inherent to the adopted T e (and to 
a much lesser extent n e ) as well as the utilized methodology. 
The most obvious example would be the map showing the dis- 
tribution of the argon abundance, where a different methodol- 
ogy was adopted to estimate the ICF depending on whether the 
[ Any 1/147 40 emission line was or was not detected (see Fig. 6 
in lPaper fl, to locate those spaxels with detection). Nevertheless, 
abundances derived in these two groups differ only by ~0. 1 dex, 
well within the expected uncertainties for this element. Also, 
12 + log(0/H) marginally anti-correlates with r e ([Om]) in the 
spaxels with lowest surface brightness. The mean (+standard de- 
viation) metallicity in our fov is 12 + \og(0/H) - 8.26 (+0.04) 
while the relative abundances \og{Ne/0) and log( Ar/O) are 
-0.64(+0.03) and -2.32(+0.05). Supporting this homogeneity, 
in all three cases - oxygen, neon and argon - the distribution of 
the values measured over the fov can be well reproduced by a 
Gaussian with <x ~ 0.03 - 0.07 dex and all the values measured 
in the individual spaxels are consistent with the mean value at 
the 3cr level. 

The abun dance in nitrogen differs from this pattern. As 
reported in iPaper il there is a roughly elliptical area of ~ 
80 pcx35 pc, associated to the main Hn region and centered 
on the two massive SSCs, that presents an enhancement of N/O. 



More interestingly, a putative secondary excess in nitrogen not 
reported so far appears in the vicinity of (but not centered on) 
knot jj3, at the area of our fov with the lowest surface brightness. 



The existence of a second area with slight N / O enhancement 
is supported by the histogram of the values measured along the 
FLAMES fov (not shown). Contrary to the case for O, Ne and 
Ar, three Gaussians centered at log(N/0) = -1.32, —1.17 and 
-0.95 and with widths u = 0.05,0.07,0.03 are needed to re- 
produce the full distribution. Most of the spaxels are associated 
to the first Gaussian and trace the typical N/O abundance for 
NGC 5253. The second largest group, with log(N/0) ~ -0.95, 
are associated to the area with N-enhancement already reported 
for the main GH nR. Using the criterion that spaxels belonging 
to two Gaussian trace different abundances if their centers are 
separated by more than 3<x, we can conclude that the third group 
of spaxels, which are associated to the area in the vicinity of 
knot #3, traces a zone with slightly larger, but clearly distinct, 
N/O value than the typical one for the galaxy. Given that this 
area is located at the spaxels with the lowest surface brightness, 
and to reject the possibility of any systematic effect due to a 
poor S/N ratio in faint lines like [Om]/14363, we extracted a 
spectrum by summing up the flux in a rectangular area of 5x6 
spaxels at this location. The measured relative N-abundance, 
log(N/0) - -1.15, is in agreement with the result found on a 
spaxel-by-spaxel basis. 
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Fig. 9. Electron density maps from different tracers. From top to 
bottom: [On], [Sn], and [Fein]. The position of the three main 
peaks of continuum emission are shown as crosses for reference. 



4. Discussion 



4.1. Electron density structure 



As we discussed in Sec. l3.3l for a given spaxel, different ions are 
associated to different layers along the line of sight. We have pre- 
sented maps for tracers of elect ron dens ity based on [S n], [Fe m] 
and [On] emission line ratios (Pa per J. and Sec. I3.2l >. The cor- 
responding density maps were derived using the same method- 
ology as in Sec. I3.1.1I and appear in Fig. [9] Values consistent 
with being below the low-density limit are found in most of the 
fov, independently of the utilized tracer, while knot #2 presents 
values somewhat larger (~190 cm- 3 ). The richest density struc- 
ture is found in the main GH nR. Here, all three tracers depict 
an n e structure with two peaks: the first one centered at knot |}1 
while the second one at ~2'/6 (~50 pc) towards the northwest. 
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Fig. 10. Electron density maps for the different components us- 
ing the results for [O n]. Top: Narrow component; Bottom: Broad 
component. The position of the three main peaks of continuum 
emission are shown as crosses for reference. 



However the values of the peaks vary depending on the utilized 
tracer. The largest densities if Fig. |9]are traced by [S n], while 
[On] and [Fern] predict densities about 10-25% smaller. Given 
their ionization potentials, (23.3, 35.1, 30.7 eV for [Sn], [On] 
and [Fein], respectively), in the absence of extinction, n c ([On]) 
- for example - traces the density i n a layer clo ser to the ion- 
izing source than n c ([Sn]) (see e.g. lGarnettlll992l) . This can be 
reversed in conditions of heavy extinction. Even if the intrinsic 
structure would be the same, lines involved in the determination 
of « e ([On]) are bluer (and therefore more sensitive to extinc- 
tion) than those associated to n e ([S n]) which therefore can probe 
deeper in the nebula. Since this GH nR suffers from relatively 
high extinction (see e.g. Fig 3. in iPaper II) the differences found 
between the three density maps are consistent with an onion- 
like structure where inner layers are denser than the outer ones. 
This is supported by the relative densities found between the two 
peaks. Also it is consistent with the fact that the [Fe m] ratio doe s 
not trace very high densities as the [Ariv] ratio did (Pa per Hi . 
With an ionization potential of 59.8 eV, [Ariv] would sample 
the densest innermost layer of this onion-like structure. 

The one Gaussian fitting approach utilized up to now is use- 
ful to have a picture of the overall density structure in the galaxy. 
However, it is known that emission lines in NGC 5253 present 
complex profiles with asymmetries tracing different kinematics 
components. For the GHnR, we tentatively measured the elec- 
tron density for what we called the "bro ad" and "narrow" com- 
ponents on a spaxel-by-spaxel basis in IPaper il using the [S n] 
lines. Results were affected by large uncertainties, mainly due to 
the deblending procedure but also to the S/N ratio of the spec- 
tra. They suggested large and uniform densities for the "broad" 
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Fig. 11. Electron temperature maps from different tracers. Top: 
[Om]; Bottom: [Sn]. The position of the three main peaks of 
continuum emission are shown as crosses for reference. 
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Fig. 12. Comparison between the predicted r e ([Om]) and 
r e ([Sn]). Spaxels were grouped in the same three bins in 
n c ([On]) as in Fig. [7] as it is indicated in the upper right-hand 
corner of the diagram. 



component (n e ~ 470 cm 3 ), while the "narrow" component pre- 
sented a decrease of density from the northwest to the south- 



east. The stronger [O n] lines offer the possibility of discussing 
this difference with much reduced uncertainties. Electron den- 
sity maps from our multi-component analysis are presented in 
Fig. [10] The structure depicted in the map for broad component 
(bottom panel), with larger densities (~ 500 cm" 3 ) close to knot 
jjl and decreasing outwards, is consistent with the general pic- 
ture sketched above (i.e. larger densities closer to the ionizing 
source). The narrow component (upper panel) present a different 
structure, with the highest density (<600 cm" 3 ) at the secondary 
peak of the GH iiR at ~ 2'.'6 from k not Hi. T his can be under- 
stood within the proposed scenario in iPaper il (see their Fig. 21, 
where different elements associated with the area of the GH nR 
are shown). There, this narrow component was associated to a 
shell of previously existing quiescent gas that has been reached 
by the ionization front. Densities in the shell would be high, if 
this pre-existing gas had been piled-up by the outflow associated 
to the broad component. 



4.2. Electron temperature structure 

As it happened with n e , a comparison of T e maps derived from 
lines associated to different ions is useful to asses the T e struc- 
ture in 3D of NGC 5253. Derived maps for r e ([Om]) and 
T e {[S n]) are presented in Fig. [TTJ At present, maps of T e based 
on any tracer are still s carce and most of the t i mes focused 
on G alactic objects (e.g. iNunez-Dfaz et al ] l2012t rTsamis et al. 
l20Q8h . Indeed, to our knowledge this is the first time that a map 
for T e ([S n]) is provided and one of the few existing examples of 
T e ([Om]) maps in extragalactic astronomy. As we pointed out 
in Sec. 13.21 the overall structure is the same for r c ([Om]) and 
r e ([Sn]) maps. However, r f ([Sn]) is smaller than r f ([Om]), 
with typical T e ([S ii])/r e ([Om]) ranging from -0.8 at the GHiiR 
to ~0.6 in the areas of low surface brightness. 

This relation between both temperatures is shown on a 
spaxel-by-spaxel basis in Fig.[J_2] Only those spaxels where the 
estimated uncertainties for both temperatures are smaller than 
40% were considered. Both, the different values for r e ([Oin]) 
and T e ([S n]) and the larger differences at low surface brightness, 
can be understood within the frame of a change in the relative 
contribution along the line of sight of the warmer gas associated 
to the GH nR and a colder and more diffuse gas component that 
extends further away. The 2D information on the plane of the 
sky (i.e. the maps in Fig. [TTJ) fits also in a satisfactory manner 
with this interpretation. 

Summarizing, both the available information on the plane of 
the sky (Fig. [TTJ) and along the line of sight (Fig.[T2b are consis- 
tent with a T e structure in 3D with higher temperatures close to 
the ionizing source surrounded by a more diffuse component of 
ionized gas at lower temperatures. 



4.3. Validity of strong line methods: N/O relative abundance 

Ideally, metallicity and relative abundances for the ionized gas 
in galaxies should be derived in a direct manner. However, this 
requires the determination of T e , which depends on the detec- 
tion of faint lines like [Om]/l4363. In extragalactic astronomy, 
most of the time this is not feasible. Instead, certain combi- 
nations of strong emission lines with more or less well estab- 
lishe d empirical and/or theoretical calibrations should be used 
(e.g. iPerez-Montero & Continill2009t iKewlev & Eliisonl 12008). 
A comparison between both methodologies on a spaxel-by- 
spaxel basis would be useful to test the reliability of the strong 
line-based methods and identify the cause of possible biases. For 
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Fig. 13. Predicted N/O maps using the N202 {top) and N2S2 
(bottom) tracers. Note that for an easier comparison with the pre- 
dictions of the direct methods, maps are displayed using a scale 
covering the same range as the one presented in Fig. [8] 



example, iMonreal-Ibero et al.l (l2011al) and iRelano et all <l201dh 
found that metallicity tracers are modulated by the ionization 
struct ure. The best-kno wn metallicity tracer would probably be 
R23 CP agel et all 1 19791) . Unfortunately, with a 12+log(0/tf) ~ 
8.28, NGC 5253 falls at the turn-over of the metallicity-R23 re- 
lation and therefore, this tracer is not appropriate for metallic- 
ity determinations in this galaxy. Alternatively, two widely used 
tracers are Q3N2 = log(([Q ml /l5007/H/?)/([ N nU6584/Hg)) 
dAlloin et all [l979h and N2=[Nn],16584/Har dDenicolo et all 
2002). However, both involve [Nn] emission lines, and are 
therefore, affected by the variation of relative N/O abundance 
across the galaxy. Thus, NGC 5253 is not an adequate exam- 
ple to test the reliability of metallicity determinations based 
on strong optical lines. Instead, since this galaxy presents a 
range of N/ O abundances, a different test for the correspond- 
ing strong line based tracers would be of high interest. This 
will be discussed in this section. As a baseline, we converted 
our line ratio maps i nto relative abundances u s ing th e ex- 
pressions provided by IPerez-Montero & Continil (|2009), that 
relate N/O with N202 = log([Nn]^6584/[On]AA3726,3729) 
and N2S2 =log([NnU6 584 /rSn1^6717,6731) as defined by 
iKewlev & Dop ita ( 2002 ) and S ab badin et all dl977h . 

Maps for the estimated relative N/O abundances are pre- 
sented in Fig. [13] They show that both tracers, N202 and N2S2, 
detect the main area with N-enhancement. Interestingly enough, 
only N202 is sensitive to the newly discovered area, associated 
to knot $3. 

However, both tracers fail to predict the correct N/O abun- 
dance. This is not unexpected since the utilized relations are 
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Fig. 14. Comparison between the residuals for the determination 
of N/O based on the two strong line methods under considera- 
tion with respect to the direct one. The locus of equal residuals is 
indicated with a black line. An orange square marks the position 
where the three methods agree. Data were divided in three bins 
of n e as in Figs. 171 and [T2l 



valid to interpret global tendencies in a statistically significant 
sample of galaxies, while N/ O abundances for an individual ob- 
ject can depart about ±0.3 and +0.5 dex from this relation for the 
N2Q2 and N2S2 tracers respect ively (see e.g. Figs. 10 and 1 1 in 
IPerez-Montero & Continill2009l) . 

Figure [14] shows a comparison of the residuals between the 
strong line methods and the direct one. Only those spaxels with 
an estimated uncertainty lower than 0.25 dex in both residuals 
were included in the comparison. The good correlation between 
the residuals is consistent with both methods being affected by 
the same factors. Also, there is an increase of the residual with 
n e . This also implies an increase with the other physical parame- 
ters since a comparison of the different maps presented through- 
out the paper shows that in general, larger n e corresponds to 
brighter areas and with larger T e and ionization strength. Given 
the spatial coincidence of these variations, disentangling the rel- 
ative role of a given physical quantity as the cause of the vari- 
ations in the residuals is not straightforward for this galaxy and 
will not be addressed here. Mean (±standard deviation) for the 
residuals are ~ 0.11(+0.03) and ~ 0.02(+0.12) for N202 and 
N2S2, respectively. Therefore, a method based on N202 is ~4 
times less sensitive to any variation of the physical/chemical 
properties than the one based on N2S2. 

In summary, our comparison of the N / O abundances derived 
using strong line methods with respect to direct measurements 
in the search for chemical inhomogeneities within a galaxy, sup- 
ports the use of that based on N202 over that based on N2S2. 
The first method is sensitive to a wider range of N/O abun- 
dances and is more stable against variations of physical condi- 
tions within the area of interest. 

4.4. On the relationship between extra N and WR stars 

There are several works in the literature using long-slit ob- 
servations that suggest a connection between nitrogen over- 
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abundance and the presence of WR stars based on the si- 
multaneous detection of WR features and a higher than the 
expected N/O in a given object (e.g . [Walsh & Roy 1987; 
Thua n et alJ I l996t iGuseva et all I2000T iPusfilnik et all 12 004: 
Izotov et ail l2006bt lHagele et alj 120061: FPerez-Montero et alj 
2oToT I Lopez-Sanchez et al.1 120121) . Also, local associations be- 
tween WR emission and N-enhancements have been found 
by means of IFS based observations (e.g. I James etail [2009) 
supporting this connection. However, the general picture is 
far from this simp l e one -to-one association. For example, 
iPerez-Montero et ail d201 lb studied a sample of BCD galaxies 
that are N-overabundant over large areas supporting the idea 
of enrichme nt caused by accreti on of less chemically processed 
gas (see also l Kehrig et al.ll2008h . A similar result was found by 
Amorin e t al.1 d2012l) for a set of BCD galaxies at higher red- 
shift. On top of that, resolution effects also play a important role. 
Already at distances of ~25-40 Mpc the whole area covered here 
would have a projected size of the order of the resolution of typ- 
ical ground based obs ervation s (~1'.'0). 

We reported in iPaper II areas with both extra-N and 
WR emission in NGC 5253, as for exam ple, our knot 
fll (w id ely discussed in t he lite rat ure, e.g. Walsh & Royl 

119971: ISchaerer et aTT ~ 
together with areas with 



19891: 



1999; 



Kobulnickv et al 
Lopez-Sanchez et al I 120071) together with areas with WR 
emission but without any chemical anomaly (e.g. our knot 
(f2). Moreover, we report in Sec. 13.41 a new area presenting an 
overabundance of nitrogen without any WR emission associated 
to it. This variety of options (i.e. WR emission and enhanced 
N, no WR emission and enhanced N, WR emission and no 
enhanced N, no WR emission and no enhanced N - i.e. the 
normal case) seems a more natural situation than a one-to-one 
association since the chemical enrichment of the warm ionized 
medium due to the material expelled by the massive stars is a 
complex process. In a simplified manner, this can be divided in 
three basic steps. 

Firstly, stars must bring their processed material to the sur- 
face. Recent stellar evolution models show how the inclusion of 
rotation favors the appearance of processe d material at stages 
as early as the main sequence (i.e. [M aeder & Meynet 2000; 
iMevnet & Maederi 120051: iMaeder & Mevnetl l2010h . However 
they are particularly p owerful at the WR , or even Luminous Blue 
Variable, phases (see ICrowtherl 120071 and references therein). 
Also, observational evid ence of N-enhancement in WR ring neb- 
ulae has been reported (Fernandez-Martin et al. 2 012b . 

Secondly, this material should be expelled via stellar winds, 
first, and supernova explosions, afterwards. In the case of star 
clusters, where stars with a variety of initial masses have a va- 
riety of evolutionary paths, the yield of the different elements 
varies in a non-trivial manner with time (e.g. lMolla & Terlevichl 
l20l2h . An extra caveat arises for clusters with stellar masses 
of <10 4 M since in this regime the Initial Mass Function is 
not properly sampled and stochastic effects are important (e.g. 
IVillaverde et aDl2D10h . 

Finally, this new material must be disseminated and then di- 
luted in the warm ionized medium which, in due time, would 
reach a new chemical homogeneity. This is an even less triv- 
ial step since on top of the specific characteristics of a given 
star cluster, factors like the existence and characteristics of 
neighbouring clusters, presence of cloudlets of gas, degree of 
the inhomogeneities, etc. are important to properly trace the 
evolution of the yielded m aterial (e.g. iTenorio-Taglel [l996; 
Ide Avillez & Mac Lowll2002l) . 

Therefore, a unique evolutionary scenario to describe the 
path from the creation of the new nitrogen to its incorporation 



to the warm ionized medium seems unlikely. In the following 
we propose some evolutionary paths for our knots $1-3. These 
scenarios do not intend to be more than reasonable suggestions 
that compile the constraints derived from our results and which 
could be tested with detailed modeling. 

Regarding our knots (J2 and |}3, we hypothesize here about 
the possibility that they constitut e two snapshots of the same 
evolutionary path. According to Har ris et al.l (12004b the main 
clusters associated to them are relatively similar in terms 
of mass and youth, although those in knot |}3 are slightly 
oldefl Supporting this youth, candidat es to supernova remnant s 
have been found close to both knots dLabrie & Pritchetl 12006). 
Moreover, only knot f|2 presented spectral features typical of 
Wolf-Rayet stars. In this evolutionary path, stars would expel 
their processed nitrogen at very early stages (<5 Myr). The dif- 
ference between the age of the cluster in knot |}3 (~10 Myr) and 
the moment when the yielding of extra -N took place (as soon as 
~ 2.5 Myr, according to the models of Moll a & Terlevichll2012l) 
sets an upper limit of ~8 Myr for the duration of this process. 
If rotating O stars were playing a relevant role, this limit could 
be extended up to the age of clusters. A lower limit can be esti- 
mated under the most optimistic (and efficient) assumption: in- 
stantaneous incorporation of the yielded material to the warm 
ionized gas. If, in addition, we assume that the veloci ty of the 
ionized gas with respect to knot |}3 (Av~20 km s _1 , see iPaperl) 
traces the velocity under which the contamination of extra ni- 
trogen propagates through the ISM, the area could be enriched 
in only ~2 Myr. Therefore, the process of the cooling down and 
mixing of the yielded material with the ISM over an area of 40- 
50 pc in diameter should last between ~2 and ~8 Myr. Note 
that the moment when the newly created material has completely 
mixed with the previously existing gas, and any chemical inho- 
mogeneity has been erased, should occur much later and cannot 
be delimited with these observations further than ~10 Myr. 

Knot )il should follow a different evolutionary path since 
at this location both WR features and nitrogen enrichment are 
found. The knot is associa ted to two very young and mas- 
sive super star clusters (e.g.lSchaereretal.lll999i iGoriian et al.l 
2001; lAlonso-H errero et al. 2004), embedded in a very dense 
and compact nebula dTurner et aTl 20001) . In that sense, they can 
be seen as nascent clusters that have not managed yet to dis- 
perse the cloud of gas where they were born. This implies a 
very particular set of physical conditions for the warm inter- 
stellar medium. Indeed, a gi ant molecular clou d associated to 
this region has been reported (Meier et al. 2002) and we showed 
through this work that n e and T e associated to this region are 
relatively high, as well as the degree of ionization. Also, the 
region p resents supersonic velocity widths and high extinction 
dPaper J) . Somehow, the combination of some of these particu- 
lar conditions has made possible the incorporation of the newly 
created nitrogen at an earlier stage. 

A third example mentio ned in previous work provides 
a diffe rent evolutionary path. iGarcfa-B enito & Perez-Monterol 
(2012) gives an age of 3.5 Myr and a stellar mass of ~ 5 x 10 3 M G 
for their knot C in NGC 6789. Moreover, they report a relative 
abundance of nitrogen ~0.2 dex larger than in their other aper- 
tures. No WR feature was detected for this knot (Garcfa-Benito, 
private communication). Clearly, at this mass regime, the yield 
is dominated by stochastic effects. Given the lack of any WR de- 



7 Knot ft2 corresponds to clusters 4 and 8 in iHarris et"aLl (120041) . 
with stellar masses of ~2.7 and 1.3xl0 4 M Q and ages of 1 and 5 Myr 
while knot (13 corresponds to their cluster 3 and 5, with masses 4.2 and 
2.1xl0 4 M G and ages of 8 and 11 Myr, respectively. 
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tection one has to resort to rotating O stars as the cause of the 
N-enhancement. Contrary to what happened with our knot (jl no 
particularly extreme physical conditions other than a complex 
inner dust structure were reported for its surrounding ISM. This 
leaves the open question of how this ISM managed to incorpo- 
rate this new material in such a short time scale. 

All in all, even if there seems to be a connection between WR 
emission and nitrogen enhancement, as supported by the fact that 
WR galaxies show an elevate d N/O relative to non-WR galax- 
ies dBrinchmann et al.l 12008). local examples like NGC 5253 
and NGC 6789, where linear spatial resolutions of ~20 pc 
arcsec -1 can be achieved, show how this relationship is complex. 
Specifically, the case of the knots $2 and f)3 in NGC 5253, where 
WR emission and N-enhancement are associated to different star 
clusters separated by only ~90 pc, illustrates the possibility that 
in galaxies at distances >25 Mpc, the spatial coincidence of WR 
emission and N/ O overabundances does not necessarily imply 
an intimate association between them based on cause-effect re- 
lationship. However, this does not reject the possibility that both 
of them were related due to a common external cause. As an ex- 
ample, in NGC 5253 one could identify as this external cause the 
putative event that triggered the starburst (i.e. the past interaction 
with M 83). 

5. Conclusions 

We have carried out a detailed 2D study of the physical and 
chemical properties of the ionized gas in the central part of 
NGC 5253, a very nearby BCD. The area was mapped in a con- 
tinuous manner with the ARGUS -IFU unit of FLAMES. This 
work re presents the natural continuation of the one presented in 
Paper I. The different maps utilized along the paper as well as the 
reduced data cubes are available as FITS files from the authors. 
The major conclusions can be summarized as follows: 

1 . Physical and chemical properties of the ionized gas associ- 
ated to the main star clusters were derived by extracting spectra 
in apertures of 9-10 spaxels. Measurements associated to knots 
(jl and ((2 agree in general with those previously reported. The 
existing discrepancies are associated to the assumed electron 
temperature. Abundances for knot (J3 (not reported so far) are 
also provided. With the exception of the relative abundance for 
N/ O in knot fll and possible knot $3, no chemical species seems 
overabundant in any of the selected apertures. 

2. Maps of the electron density based on four different trac- 
ers - namely [On], [Sn], [Fern], and [Anv] line ratios - were 
discussed. In all the cases, higher densities are associated to the 
main GHnR. The joint analysis of these maps is consistent with 
a 3D stratified view of the nebula where the highest densities are 
located in the innermost layers while density decreases when go- 
ing outwards. 

3. The 2D n e structure for the two main kinematic compo- 
nents was also derived using the [O n] lines as baseline. While 
the so-called broad component follows the picture described 
above, there is a change of structure for the n arrow one. This 
fits well with the proposed scenario in Paper I where this com- 
ponent was associated to a shell of previously existing material 
that had been piled up by the outflow associated to the broad 
component. 

4. Maps for ^([Sn]) and ^([Om]) were derived. To our 
knowledge, this is the first time that a T e {[S n]) map for an ex- 
tragalactic object is presented. Also, we provided with one of 
the few examples of existing r e ([Om]) maps up to date. The 
joint interpretation of the information on the plane of the sky 
and along the line of sight is consistent with a T e structure in 



3D with higher temperatures close to the main ionizing source 
surrounded by a colder and more diffuse component. This T e 
structure together with the lack of any strong broad component 
far from the main GH nR dPaper J) is in accord with the lack of 
any clear evidence of shocks playing a dominant role. 

5. Ionization structure was traced by means of the 
[Ora]Xl4959,5007/[On]/l/l3726,3729, [On]^3726,3729/H^, 
and [S n]/M6717,6731/Har ratios. The two first of them predict 
similar ionization degree while the third one would be typical 
of lower ionization. A possible 3D interpretation of both the 
observed structure in the maps for each individual ratio and the 
discrepancy of [S n]AA67 17,673 1/Ha for individual spaxels is 
consistent with the lower ionization species (i.e. S + ) delineating 
the more extended diffuse component. 

6. Maps for the 2D distribution of abundances for oxygen, 
neon, argon were derived. All of them are consistent with no 
chemical inhomogeneities. The derived mean (+ standard devi- 
ation) oxygen abundance is 12 + \og(0/H) = 8.26 + 0.04. The 
mean logarithmic relative abundances for argon and neon were 
-0.65 ± 0.03 and -2.33 + 0.06, respectively. 

7. In the same manner, a map for the 2D distribution of ni- 
trogen was derived. NGC 5253 typically presents a log(N/0) 
of ~ -1.32 + 0.05. However, there are two locations with en- 
hanced N/ O. The first one was already reported and character- 
ized. With a log{N/0) ~ -0.95, it occupies an elliptical area 
of about 80 pcx35 pc and is associated to the two SSCs at the 
nucleus of the galaxy. The second one is reported here for the 
first time. It presents a log(N/0) ~ —1.17 and it is associated 
to two moderately massive (2 - 4 x 10 4 M ) and relatively old 
(~ 10 Myr) clusters (knot (j 3). 

8. The map of N/O relative abundance derived through the 
direct method was compared with those derived using strong line 
methods. The comparison supports a method based on N202 
over a method based on N2S2 in the search of chemical inho- 
mogeneities within a galaxy since the first method is sensitive 
to a wider range of N/ O abundances and is more stable against 
variations of physical conditions within the area of interest. 

9. We utilized the results on the localized detection of WR 
emission and nitrogen enhancement to compile and discuss the 
factors that affect the complex relationship between the presence 
of WR stars and N/O excess. Even if there seems to be such a 
relationship, a unique scenario describing the path from the pro- 
duction of the new nitrogen to its incorporation into the warm 
ionized medium seems unlikely. In particular, we use the areas 
associated to knots $2 and (j3 in NGC 5253 as examples of WR 
emission and N-enhancement that would be perceived as spa- 
tially coincident at distances >25 Mpc but that are not intimately 
associated in a cause-effect fashion. However, this does not re- 
ject the possibility that both of them are related due to a common 
external cause. 
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